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The algorithm assumes a constant bit rate over a timing window of specified duration (e.g., a 
specified number of consecutive frames), where the current frame to be encoded lies in the interior of 
the timing window. A target bit rate for the current frame is initially selected by calculating the number 
of bits already used to encode other frames within the window and then assuming that the remaining 
available bits allocated to the timing window will be evenly distributed to the remaining unencoded 
frames in the timing window. The target bit rate may then be optionally adjusted based on scene 
content, encoder state, and buffer considerations. Through a combination of target bit allocation and 
frame skipping, spatial and temporal resolutions are maintained within acceptable ranges while meeting 
buffer delay constraints. The algorithm has also been extended to include PB frames in addition to 
P-only coders. 



10^ 



cucuun nnti or bits wed 

10 QCDOt MKT lit fWrtS 

ii neat 



104^ 



ULCIUU mWtt OF UTS 
AMLAJU II BEST Of 1IM0W 



106^ 



cucuun war t assuim 
equal *mn or bib/fume 
fat wsr of i iva 



108- 



A0JUST TARGET T USED OH 
QUUTIWnX SATWATK* 



ADJUST 


TUBE! T BASED Of 


CURHBT 


aid ria mnairtou 



AOJtST TARGET T BASE 


0 01 


lUBEA OF BITS USED 


FOR 


ntETIOUS MAKE 





ADJUST TIMET T BASED OH 
Biff El flXLHESS 



ADJUST TARGET T- BASED CN 

mwia ot high ton* 



118- 



14 OC0C! CUttOT FRAKE 
BASED a tAflST r 



FOR THE PURPOSES OF INFORMATION ONLY 



Codes used to identify States party to the PCT on the front pages of pamphlets publishing international applications under the PCT. 



AL 


Albania 


ES 


Spain 


LS 


Lesotho 


SI 


Slovenia 


AM 


Armenia 


FI 


Finland 


LT 


Lithuania 


SK 


Slovakia 


AT 


Austria 


FR 


France 


LU 


Luxembourg 


SN 


Senegal 


AU 


Australia 


GA 


Gabon 


LV 


Latvia 


sz 


Swaziland 


AZ 


Azerbaijan 


GB 


United Kingdom 


MC 


Monaco 


TD 


Chad 


BA 


Bosnia and Herzegovina 


GE 


Georgia 


MD 


Republic of Moldova 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagascar 


TJ 


Tajikistan 


BE 


Belgium 


GN 


Guinea 


MK 


The former Yugoslav 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Greece 




Republic of Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungary 


ML 


Mali 


TT 


Trinidad and Tobago 


BJ 


Benin 


IE 


Ireland 


MN 


Mongolia 


UA 


Ukraine 


BR 


Brazil 


IL 


Israel 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


ES 


Iceland 


MW 


Malawi 


US 


United States of America 


CA 


Canada 


IT 


Italy 


MX 


Mexico 


uz 


Uzbekistan . 


CF 


Central African Republic 


JP 


Japan 


NE 


Niger 


VN 


Viet Nam 


CG 


Congo 


KE 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


CH 


Switzerland 


KG 


Kyrgyzstan 


NO 


Norway 


ZW 


Zimbabwe 


CI 


Cote d'lvoire 


KP 


Democratic People's 


NZ 


New Zealand 






CM 


Cameroon 




Republic of Korea 


PL 


Poland 






CN 


China 


KR 


Republic of Korea 


PT 


Portugal 






CU 


Cuba 


KZ 


Kazakstan 


RO 


Romania 






CZ 


Czech Republic 


LC 


Saint Lucia 


RU 


Russian Federation 






DE 


Germany 


LI 


Liechtenstein 


SD 


Sudan 






OK 


Denmark 


LK 


Sri Lanka 


SE 


Sweden 






EE 


Estonia 


LR 


Liberia 


SC 


Singapore 







WO 00/18137 



PCT/US99/21841 



FRAME-LEVEL RATE CONTROL FOR VIDEO COMPRESSION 
BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to image processing, and, in particular, to video compression. 

Cross-Reference to Related Applications 

This application claims the benefit of the filing date of U.S. provisional application no. 
60/100,939, filed on 09/18/98 as attorney docket no. SAR 12728P. 

Description of the Related Art 

The goal of video compression processing is to encode image data to reduce the number of bits 
used to represent a sequence of video images while maintaining an acceptable level of quality in the 
decoded video sequence. This goal is particularly important in certain applications, such as real-time 
video conferencing, where transmission bandwidth limitations may require careful control over the bit 
rate, that is, the number of bits used to encode each image in the video sequence. In order to satisfy the 
transmission and other processing requirements of a video conferencing system, it is often desirable to 
.have a relatively steady flow of bits in the encoded video bitstream. 

Achieving a relatively uniform bit rate can be very difficult, especially for video compression 
algorithms that encode different images within a video sequence using different compression 
techniques. Depending on the video compression algorithm, images may be designated as the 
following different types of frames for compression processing: 

o An intra (I) frame which is encoded using only intra-frame compression techniques, 
o A predicted (P) frame which is encoded using inter-frame compression techniques based on a 
previous I or P frame, and which can itself be used as a reference frame to encode one or more other 
frames, 

o A bi-directional (B) frame which is encoded using bi-directional inter-frame compression 
techniques based on a previous I or P frame and a subsequent I or P frame, and which cannot be used 
to encode another frame, and 

o A PB frame which corresponds to two images - a P frame and a subsequent B frame - that are 
encoded as a single frame as in the H.263 video compression algorithm. 

Depending on the actual image data to be encoded, these different types of frames typically require 
different number of bits to encode. For example, I frames generally require the greatest numbers of 
bits, while B frames generally require the least number of bits. 

In a typical transform-based video compression algorithm, a block-based transform, such as a 
discrete cosine transform (DCT), is applied to blocks of image data corresponding to either pixel 
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values or pixel differences generated, for example, based on a motion-compensated inter-frame 
differencing algorithm. The resulting transform coefficients for each block are then quantized for 
subsequent encoding (e.g., run-length encoding followed by variable-length encoding). The degree to 
which the transform coefficients are quantized (also referred to as the quantization level) directly 
5 affects both the number of bits used to represent the image data and the quality of the resulting decoded 
image. In general, higher quantization levels imply fewer bits and lower quality. As such, quantization 
level is often used as the primary variable for controlling the tradeoff between bit rate and image 
quality. 

At times, using quantization level alone may be insufficient to meet the bandwidth and quality 
10 requirements of a particular application. In such circumstances, it may become necessary to employ 
more drastic techniques, such as frame skipping, in which one or more frames are dropped from the 
video sequence. Such frame skipping may be used to sacrifice short-term temporal quality in the 
decoded video stream in order to maintain a longer-term spatial quality at an acceptable level. 

15 SUMMARY OF THE INVENTION 

The present invention is directed to an algorithm for achieving frame-level bit-rate control for 
video compression processing for applications such as video conferencing in which real-time 
processing and transmission requirements make such rate control especially important. According to 
one embodiment, the present invention is a method for controlling bit rate by a machine implementing 

20 video compression processing of a video sequence to select a target bit rate for a current frame in the 
video sequence, comprising the steps of (a) generating a number of bits used to encode a previous M 
frames in the video sequence, wherein M is a specified number of frames; (b) generating a number of 
available bits for an unused portion of a multi-frame bit budget corresponding to a specified number of 
bits available for encoding N frames in the video sequence, wherein N is a specified number of frames 

25 and M < N, by subtracting the number of bits used to encode the previous M frames from the multi- 
frame bit budget; (c) generating the target bit rate for the current frame by dividing the unused portion 
of the multi-frame bit budget by (N-M); and (d) encoding the current frame based on the target bit rate. 

BRIEF DESCRIPTION OF THE DRAWINGS 
30 Other aspects, features, and advantages of the present invention will become more fully 

apparent from the following detailed description, the appended claims, and the accompanying drawings 
in which: 

Fig. 1 shows a flow diagram of the processing implemented by a P-frame-only coder to select a 
target number of bits for encoding the current frame, according to one embodiment of the present 
35 invention; and 
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Fig. 2 shows a flow diagram of the processing corresponding to a frame-skipping strategy, 
according to one embodiment of the present invention. 



DETAILED DESCRIPTION 
In one embodiment, the present invention is directed to a frame-level rate control and skipping 
strategy for a real-time video conferencing system that allocates frame-level targets (i.e., allocations of 
bits) based on scene content, encoder state, and buffer considerations. The algorithm maintains a 
constant bit rate over a timing window of specified duration (e.g., specified number of consecutive 
frames), while allowing variations in the target and actual bit rates used to encode individual frames. 
Through a combination of target bit allocation and frame skipping, spatial and temporal resolutions are 
maintained within acceptable ranges while meeting buffer delay constraints. The algorithm may be 
extended to include PB frames in addition to P frames. - 

Notation 

The following notation is used in this specification: 
R: Channel transmission rate in bits/sec (e.g., 24000 bits/sec). 

B PP : Bits per picture = R/original_frame_rate (e.g., for an original frame rate of 30 

fps and R=24Kbps, Bpp=800 bits per picture). 

Bprev: Bits used to encode the previous frame. 

fs: Frame skip corresponding to the desired average frame rate (e.g., for input 

video having 60 fps and a transmission rate of 15 fps, fs=4). 

act_buf_size: The actual size of the encoder buffer. This is governed by the picture format. 

VB Vf_b: Buffer fullness before putting the currently encoded frame bits into buffer. 

VBVf_a: Buffer fullness after putting the currently encoded frame bits into buffer. 

max_buf_size: The maximum size of the steady state buffer. In one embodiment, this is 
chosen to be R/2, which means that 0.5 seconds is the upper bound on buffer delay. 

ssbs: The desired VBVf_a (steady state buffer size). 

S, Sp: Interframe distortion after motion compensation for the current frame and 

average distortion over a specified number of previous frames, respectively. 

Frame Skipping Strategy 

In order to meet buffer delay constraints while maintaining the spatial quality within acceptable 
limits at a given bit rate, frame skipping may be adopted. In one embodiment, the average delay for a 
frame in the buffer is constrained to be ssbs/R, and the maximum delay is dictated by the accuracy with 
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which the target is met. Usually, it is only marginally larger that ssbs/R. The skip after encoding a 
picture is computed as follows: 

Set skip=l , and VB VfJ> = VB Vf_a - Bpp. 
while (VBVf_b + Bf > ssbs){ 
5 if (VB Vf_b - Bpp < 0) break; /* this avoids underflow and under-utilization of bits */ 

skip++ 

VB Vf_b = VBVf_b - Bpp 

} 

Bf is the estimate for the next frame to be encoded and is computed as (Bprev + Bpp*fs)/2. Thus, the 
10 algorithm attempts to maintain the buffer delay around ssbs/R. 

Frame-Level Rate Control for P-Frame-Onlv Coders 

Fig. 1 shows a flow diagram of the processing implemented by a P-frame-only coder to select a 
target number of bits (T) for encoding the current frame, according to one embodiment of the present 
15 invention. The processing depicted in Fig. 1 is described in the following sections. 

Target Computation over a Window of Frames 

Under the processing of Fig. 1, an attempt is made to maintain a constant bit rate over a 
window of W frames, centered at the current frame being encoded. For the first Nl frames in the 
20 window, the number of bits (B 1) used to encode those frames is calculated (step 102 of Fig. 1). The 
number of frames left to be encoded is N2=(W-Nl)/fs, where Nl is the total number of frames 
processed including those that are skipped. The available number of bits for the whole window is Bw 
= W*Bpp. The target T for the current frame is T = (Bw-Bl)/N2 (steps 104 and 106). This 
computation comes from the assumption that the same number of bits will be used for each of the 
25 remaining frames in the window. 



Target Modification 

This section describes an algorithm to modify the target T based on the state of the encoder and 
the scene content. The target is modified based on one or more of the following optional 
30 considerations: (a) quantizer saturation, (b) distortion after motion compensation as compared to 
distortion over the previous frames, (c) bit count for the previous frame, and (d) buffer fullness. 

In particular, for quantizer saturation, if the quantizer parameter QP saturated for the previous 
picture (e.g., mean QP >25, for a system having 31 quantization levels), the current frame target T may 
be increased (step 108 of Fig. 1) as follows: 
35 T=T*(1 + p*(QP - QP.thresh)) 
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where (3 is a specified factor (e.g., 0.06) and QP_threshold is a specified saturated quantization level 
(e.g. 25). 

While the quantizer parameter is not saturated, the target T may be modified based ori current 
and past distortions (step 110) as follows: 

T=T*(Sp + k*S)/(k*Sp + S), 
where S is the motion-compensated distortion (e.g., basedon a sum of absolute pixel differences 
(SAD) measure) for the current frame, Sp is the average motion-compensated distortion over the past 
M frames (e.g., M = 2), and k>l (e.g., k=4). Sp is preferably a weighted average that emphasizes more 
recent data. 

The target T may also be modified based on the number of bits used for the previous frame 
(Bprev) (step 112) as follows: 

T=a*T+(l-oc)*Bprev, 

where 0<=a<=l (e.g., a =0.7). This step reduces fluctuation in the quality of reference frames, and 
helps ensure that any degradation over time is graceful. 

The target T may also be modified based on buffer fullness (step 114) as follows: 

T=T*(a+A,*b)/(A,*a+b), 

where a = VBVfb (the buffer fullness before putting the currently encoded frame bits into the buffer), 
b = ssbs - VBVf_b, and A>1. This enables the algorithm to take advantage of the buffer. If VBVf_b is 
less than half of the desired steady state buffer size (i.e., ssbs/2), then the target is increased; otherwise, 
the target is decreased. 

The resulting target T is then used to encode the frame (step 118), for example, using 
macroblock adaptation of the QP (quantization parameter) based on the TMN 5 orTMN 8 strategies 
for the H.263+ standard. 

High-Motion Exception 

In high-motion regions, the above algorithm may lead to loss of spatial quality and the 
algorithm might take a long time to recover. Thus, in these situations, the algorithm may take 
advantage of the buffer and choose to boost the target (step 116 in Fig. 1). Furthermore, when the 
algorithm returns to "steady-state " the algorithm gracefully decreases the frame target to the steady 
state value (i.e., approximately fs*Bpp). Graceful degradation may be achieved based on a weighted 
average on the target with the previous bit count Bprev, for example: 

T = 0.8T + 0.2Bprev 

This would enable the algorithm to maintain the spatial quality of the frames even in the presence of 
abrupt motion, and allows the algorithm to recover faster when the scene motion reduces. 
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The high-motion detection may be based on a comparison between the current motion- 
compensated distortion (S) and the weighted average of the previous motion-compensated distortions 
(Sp). In one embodiment, the target modification is implemented as three linear segments selected to 
approximate a smooth curve, where: 
5 T = (1 + 4*factor)/(4 + factor)*T for factor < 1.05 

T = (1.03 + (factor-1.05)/0.3)*T for 1.05 < factor < 1.15 
T = (1.36 + 4*(factor-1.15))*T for factor > 1.15 
where factor = S/Sp. 

10 Frame-Level Rate Control for Coders using PB or P Frames 

The strategy described above was designed to code P frames only. In this section, the strategy 
is extended to coders with PB frames. 

When P and PB frames are mixed, the delays for the P frames are significantly lower than the 
delay for the B frame of a PB frame (for the same skip). For ease of discussion, it is assumed that the 
15 algorithm operates at a constant frame skip fs. In that case, the overall system delay for a typical P- 
only frame from the moment of capture to availability at the decoder is given as follows: 

fs/ref_frame_rate + VBVf_b/R + processing delay + channel delay, 
while the delay for the B frame of a typical PB frame is given as follows: 

3fs/ref Jrame_rate + VBVf_b/R + processing delay + channel delay 
20 In general, VB Vf_b is lower for the PB frame than the P frame, since a PB frame has to wait for the P 
frame (which is typically 2*fs from the previous reference frame and thus depleting the buffer more) 
before encoding. 

Thus, when P and PB frames are mixed, delays for the P frames are smaller than those for the 
PB frames. To keep the delays equal, the following guidelines may be adopted: 
25 1 . The VBVf_b before encoding a P-only frame is fs*Bpp more than the VBVf_b before 

encoding a PB frame. 

2. Each P frame is delayed by fs*Bpp/R at the decoder before it is displayed. 
Fig. 2 shows a flow diagram of the processing corresponding to a frame-skipping strategy 
based on the above guidelines, according to one embodiment of the present invention. According to 
30 this strategy: 

1 . First, determine the skip after the previous reference frame based on the strategy used in 
the P-only case (step 202 of Fig. 2). The only difference comes in the case where the previous frame 
was a PB frame. In that case, the number of bits used to encode the previous frame (Bprev) is twice 
the number of bits required for one frame, and this is divided by 2 when the number of bits for the 
35 current frame is estimated. 
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2. Then, perform a PB decision to decide whether or not to code the current frame as a P 
frame. The strategy is based on the observation that H.263 PB frames do not work well for large 
motion. (This statement is not necessarily true in the improved PB mode). Thus, if large motion is 
detected (steps 204-208), the decision is made to encode the current frame as a P frame (step 210). 
One possible high-motion detection algorithm is described in the next section. If the algorithm chooses 
to encode this frame as a P frame, then the processing returns to step 1 after encoding the P frame. 

3 . If a PB frame is coded, the algorithm determines the skip between the B frame and the P 
frame of the PB frame based on the buffer fullness (step 212). Essentially, the algorithm determines 
the skip as follows: 

VBVJullness + max(Bprev, 2*fs*Bpp) < ssbs 
According to this equation, the buffer fullness after encoding the PB-frame should be less than the 
steady-state-buffer size ssbs. The term max(Bprev, 2*fs*Bpp) is used to estimate the number of bits 
for the frame. Thus, skip is incremented until the buffer fullness, after encoding, drops below ssbs. 
Then, encode the PB frame (step 214) and go to step 1 . 

In the case where high motion is detected between the B and P frames of a PB frame, the PB 
mode is turned off for this frame and the P frame of the "PB" frame is encoded as a P-only frame. This 
situation should not occur very often in a sequence. 

Motion Detection for the PB decision 

The PB decision described previously requires a high-motion detector. This detector may be 
implemented as follows. 

1 . Determine the number of blocks that are moving (Ml ) (step 204 of Fig. 2). This may be 
determined from the raw frame differences or from the motion-compensated frame differences and the 
motion field. If motion estimation is performed, Ml refers to the number of blocks having a non-zero 
motion vector. Otherwise, Ml refers to. the number of blocks having a distortion measure greater than 
some specified threshold level (e.g., SAD>100). 

2. Determine the number of blocks that have high motion (M2) (step 206), again based on 
either the raw frame differences or the motion-compensated frame differences and the motion field. If 
motion estimation is performed, M2 refers to the number of blocks whose motion vector has a 
magnitude greater than a specified threshold level. If motion estimation is not used, M2 refers to the 
number of blocks having a distortion measure greater than some specified high-motion threshold level 
(e.g., SAD>750). 

3 . If M2/M 1 is greater than a specified threshold level (e.g. , 1 . 1 5) (step 208), then code the 
current frame as a P frame (step 210); otherwise, code the current frame as a PB frame (steps 212-214). 



Summary and Conclusions 
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The frame-level rate control of the present invention is particularly applicable in real-time very 
low bit-rate coders. The algorithm can be implemented by coders using only P frames as well as 
coders using PB-frames, an added functionality when compared to the TMN8 test model for the 
H.263+ standard. The algorithm provides a buffer delay variable, which can be selected by a user to 
5 trade delay for graceful change in spatial quality over time. By adapting the target bit allocation for a 
frame to the scene content, spatial quality is maintained in high-motion areas and quick recovery is 
made possible after an abrupt motion, while maintaining the buffer delay within desired limits. 

The present invention can be embodied in the form of methods and apparatuses for practicing 
those methods. The present invention can also be embodied in the form of program code embodied in 

10 tangible media, such as floppy diskettes, CD-ROMs, hard drives, or any other machine-readable 

storage medium, wherein, when the program code is loaded into and executed by a machine, such as a 
computer, the machine becomes an apparatus for practicing the invention. The present invention can 
also be embodied in the form of program code, for example, whether stored in a storage medium, 
loaded into and/or executed by a machine, or transmitted over some transmission medium, such as over 

15 electrical wiring or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the 
program code is loaded into and executed by a machine, such as a computer, the machine becomes an 
apparatus for practicing the invention. When implemented on a general-purpose processor, the 
program code segments combine with the processor to provide a unique device that operates 
analogously to specific logic circuits. 

20 It will be further understood that various changes in the details, materials, and arrangements of 

the parts which have been described and illustrated in order to explain the nature of this invention may 
be made by those skilled in the art without departing from the principle and scope of the invention as 
expressed in the following claims. 
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CLAIMS 

What is claimed is: 

1. A method for controlling bit rate by a machine implementing video compression processing of 
a video sequence to select a target bit rate for a current frame in the video sequence, comprising the 
steps of: 

(a) generating a number of bits used to encode a previous M frames in the video sequence, 
wherein M is a specified number of frames; 

(b) generating a number of available bits for an unused portion of a multi-frame bit budget 
corresponding to a specified number of bits available for encoding N frames in the video sequence, 
wherein N is a specified number of frames and M < N, by subtracting the number of bits used to 
encode the previous M frames from the multi-frame bit budget; 

(c) generating the target bit rate for the current frame by dividing the unused portion of the multi- 
frame bit budget by (N-M); and 

(d) encoding the current frame based on the target bit rate. 

2. The invention of claim 1 , wherein step (c) comprises the further step of adjusting the target bit 
rate based on presence of quantizer saturation. 

3. The invention of claim 2, wherein step (c) comprises the step of increasing the target bit rate if 
a quantizer parameter for a previous frame was greater than a specified threshold level. 

4. The invention of claim 1 , wherein step (c) comprises the further step of adjusting the target bit 
rate based on current and past distortion levels. 

5. The invention of claim 1 , wherein step (c) comprises the further step of adjusting the target bit 
rate based on the number of bits used to encode a previous frame. 

6. The invention of claim 1 , wherein step (c) comprises the further step of adjusting the target bit 
rate based on buffer fullness. 

7. The invention of claim 1 , wherein step (c) comprises the further step of adjusting the target bit 
rate based on a characterization of motion. 

8. The invention of claim 7, wherein step (c) comprises the step of adjusting the target bit rate as 
a function of S/Sp, where S is a parameter equal to a motion-compensated distortion level for the 
current frame and Sp is a parameter equal to an average motion-compensated distortion level over a 
previous number of frames. 

9. The invention of claim 1, wherein step (c) further comprises the steps of: 
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2 (1) adjusting the target bit rate based on presence of quantizer saturation; 

3 (2) adjusting the target bit rate based on current and past distortion levels; 

4 (3) adjusting the target bit rate based on the number of bits used to encode a previous frame; 

5 (4) adjusting the target bit rate based on buffer fullness; and 

6 (5) adjusting the target bit rate based on a characterization of motion. 

1 10. The invention of claim 9, wherein: 

2 step (c)(l ) comprises the step of increasing the target bit rate if a quantizer parameter for a previous 

3 frame was greater than a specified threshold level; 

4 step (c)(2) comprises the step of adjusting the target bit rate T as follows: 

5 T=T*(Sp + k*S)/(k*Sp + S), 

6 where S is a motion -compensated distortion level for the current frame, Sp is a weighted average 

7 motion-compensated distortion over a previous number of frames, and k is a specified parameter >1 ; 

8 step (c)(3) comprises the step of adjusting the target bit rate T as follows: 

9 T=a*T+(l-a)*Bprev, 

10 where Bprev is the number of bits used to encode the previous frame and a is a specified parameter 

11 such that 0<=oc<=l; 

12 step (c)(4) comprises the step of adjusting the target bit rate T as follows: 

13 T=T*(a+X*b)/(^*a+b), 

14 where a is a parameter equal to the buffer fullness before putting currently encoded frame bits into the 

15 buffer, b is a parameter equal to ssbs-a, where ssbs is a parameter equal to a steady state buffer size, 

16 and X is a specified parameter >1 ; and 

17 step (c)(5) comprises the step of adjusting the target bit rate as a function of S/Sp, where S is a 

18 parameter equal to a motion-compensated distortion level for the current frame and Sp is a parameter 

19 equal to an average motion-compensated distortion level over a previous number of frames, wherein 

20 the function is implemented as a plurality of linear segments relating S/Sp to a factor used to adjust the 

2 1 target bit rate. 
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FIG. 1 
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